#! /usr/bin/env python

import pyx
import pyfeyn
from pyfeyn.user import *

fd = FeynDiagram()
processOptions()

in1 = Point(-4,  2)
in2 = Point(-4, -2)
in11 = Point(0,  2)
in22 = Point(0, -2)

out1= Point(4,2)
out2=Point(4,-2)

out = Point(4, 0)

vtx = Vertex(0, 0, mark=CIRCLE)

Gluon(in1,in11).addLabel(r"\Pgg")
Gluon(in2,in22).addLabel(r"\Pgg")
Fermion(in11,out1).addArrow().addLabel(r"\Pqt/\Pqb")
Fermion(in22,out2).addArrow().addLabel(r"\Pqt/\Pqb")
#Fermion(in11,in22).addArrow().addLabel(r"\Pqt/\Pqb")
Fermion(in11,in22)
Higgs(vtx,out).addLabel(r"\PHiggs")


fd.draw("./basoc_Higgs_production_Feynman.pdf")

